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Abstract 

We consider the repulsive Hubbard model on a class of lattices or graphs for which there is a large de- 
generacy of the single-particle ground states and where the projector onto the space of single-particle ground 
states is highly reducible. This means that one can find a basis in the space of the single-particle ground 
states such that the support of each single-particle ground state belongs to some small cluster and these clus- 
ters do not overlap. We show how such lattices can be constructed in arbitrary dimensions. We construct 
all multi-particle ground states of these models for electron numbers not larger than the number of localised 
single-particle eigenstates. We derive some of the ground state properties, esp. the residual entropy, i.e. the 
finite entropy density at zero temperature. 

1 Introduction 

The physics of strongly-correlated Fermi systems is one of the most exciting branches of condensed matter 
theory. The most extreme case of strong correlations occurs in systems where the single-particle problem has a 
large degeneracy. Such systems with flat bands have been studied intensively during the last 20 years. One of 
the more prominent examples is the Hubbard model with a lowest flat band. The first examples were discussed 
by Tasaki and by the present author 20 years ago 0][2][3]. ^ e a i m °f me authors was a rigorous, complete 
classification of the ground states. Among other results, it was shown that under certain conditions the Hubbard 
model has a unique ferromagnetic ground state. The final goal was the proof of the existence of metallic 
ferromagnetism in these and similar models. 

If a translationally-invariant system has a flat band, not only the Bloch states but also the Wannier states are 
eigenstates of the single-particle Hamiltonian. This means that there are localised single-particle eigenstates. 
Some of the physical properties of these systems can be related to the existence of localised eigenstates. Most of 
the proofs used to derive exact results for ferromagnetism in Hubbard models with flat bands rely on the fact that 
such localised basis states exist. On the other hand, the existence of metallic ferromagnetism, which was always 
one of the goals, is not related to localised eigenstates. One class of models where one could have metallic 
ferromagnetism is constructed from the flat band models with additional hopping terms to lift the degeneracy 
of the degenerate band [4J[3]. Another class of models have partially flat bands and thus the localised Wannier 
states are not eigenstates |6]|7). 

In the models mentioned so far the localised states, when they exist, overlap. To be precise, the representation 
of the projector onto the subspace of all degenerate single-particle ground states in position space, the single- 
particle density matrix, is irreducible. This was essential for the proofs in O |7] [6) . In Q it was shown that if 
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and only if the single-particle density matrix is irreducible, the multi-particle ground state at a special density, 
where each of the single-particle eigenstates is filled with exactly one electron, is ferromagnetic and unique up 
to the usual (2S + 1 ) fold spin degeneracy. 

Other models which recently caught some interest have strictly localised eigenstates in the sense that different 
localised single-particle eigenstates do not overlap, see e.g. Batista and Shastry [8 |, Dzerkho et al. Maksy- 
menko et al. [10|, some of the examples by Schmidt et al. [11], and the references therein. These authors 
consider examples of decorated lattices in one or two dimensions. In these systems, the physical properties are 
strongly influenced by the existence of localised eigenstates. The system has a high degeneracy for the multi- 
particle ground states. The entropy density at zero temperature is finite. Since the single-particle eigenstates 
are localised the system is likely to be non-metallic. In some cases the ground states are Wigner crystals. The 
system is paramagnetic, not ferromagnetic. 

The goal of the present paper is to provide a complete description of the class of lattices with degenerate single- 
particle eigenstates which fall into non-overlapping subsets. This is done using the projector onto the subspace 
of degenerate single-particle ground states in position space, the single-particle density matrix. The class of 
lattices with single-particle eigenstates which fall into non-overlapping subsets has a reducible single-particle 
density matrix, in contrast to the models discussed in |7, 6|. We discuss some of the properties of the Hubbard 
model on such lattices and we give a large class of examples of such lattices which can be constructed explicitly. 
The examples in |[8ll9l [T0ll belong to this class as well. The construction is possible in arbitrary dimensions. 
We construct explicitly all ground states of these models for sufficiently low densities of states - the flat band 
must be at most half filled - and we calculate the entropy density at zero temperature. We prove that there is no 
long-range order in these models. 

The paper is organised as follows: In Sect. |2]we define the class of lattices we are looking at. They are defined 
by some properties of the projector onto the single-particle ground states of the system. We give a general 
description how examples of such lattices can be constructed explicitely in arbitrary dimensions. In the Sect. [3] 
we state and proof our main result concerning the multi-particle ground states in such models. Sect. @]contains 
a summary and an outlook. 



2 Classification of the single-particle problem 

In this paper we consider a general fermionic Hubbard model 

H = Hhop + Hi nt (1) 

where 

Hhop = ^ txyC xa Cy a (2) 

{x,y}eE,a 

and 

#int = £ U x n^n xl (3) 
xeV 

on a lattice or, more generally, on a connected graph G = (V,E) with a set of vertices V and edges E connecting 
the vertices, txy are the hopping matrix elements, U x > is the local repulsive interaction. Two vertices x and y 
are connected by an edge e = {x,y} if and only if ^ 0. In this section, we consider first the single-particle 
problem in order to define the class of models (or lattices) we are dealing with. 

We consider the case where //hop has a highly degenerate single-particle ground state with eigenenergy e c /. The 
degeneracy is Nd- G does not need to be translationally invariant. In the case of a translationally-invariant 
lattice, we assume that the system has at least one degenerate energy band at the bottom of the spectrum. 
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Let B = {\ffi(x), i = 1 . . .N c /} be an arbitrary orthonormal basis in the subspace of the degenerate lowest eigen- 
states of the matrix T = (t xy )x,yev- We assume that t xy are real, however a generalisation to complex t xy is straight 
forward. Complex t xy have also been discussed in the context of flat bands, see e.g. 02 1. For real txy we choose 
the basis B to be real as well. The single-particle density matrix of these states is 

N d 

Pxy = Y,Yi(x)Yi(y). (4) 

(=1 

P = (Pxy)x yeV is the projector onto the space spanned by the single-particle ground states in position space. We 
showed that if p = (p™)* yeV is irreducible, the Hamiltonian has ferromagnetic multi-particle ground states and 
that for special particle numbers N e — Nd, the ferromagnetic ground state is unique up to the degeneracy due to 
the SI/ (2) spin symmetry Q. 

In this paper, we consider the case where the single-particle density matrix p is highly reducible, p should have 
the following properties: 

1. p is reducible. It can be decomposed into N r irreducible blocks p^, k = l,...,N r . N r should be an 
extensive quantity, i.e. N r °e °c |V|, so that in the thermodynamic limit the density of degenerate 
single-particle ground states and the density of irreducible blocks are both finite. 

2. Let Vk be the support of p,t, i.e. the set of vertices for which at least one element of p^ does not vanish. 
Pkjcy = if jc ^ Vj or y ^ Vfc. One has Vjt n Vjf = if k / k' because of the fact that p,t are irreducible blocks 
of the reducible matrix p and Uft^* Q V. 

3. We choose the B such that the support of each basis states y//(jr) is a subset of exactly one Vj.. We denote 
the number of states belonging to the cluster as v^. One has V* = Nd- 

4. V max = rnax<.{v,t} is 0(1), i.e. not an extensive quantity. 

If G represents a translationally-invariant lattice, only one or a few blocks belong to one elementary cell and the 
Vk belong to classes where, within one class, all are the same due to translational invariance. 

On each block, since p^. is irreducible, the results obtained in |7| apply. 

2.1 Lattices with such properties 

The lattices in [8, 9, 10] have the properties mentioned above. We now give a more general construction for a 
large class of lattices in arbitrary dimensions which have these properties. Let us mention that these are only 
examples and that many other lattices with reducible p^ exist. 

Our starting point for the construction of a large class of lattices or graphs with these properties is an arbitrary 
lattice or graph G — (V ,E). V is the set of vertices of G. E is the set of edges of G. We consider only simple 
graphs, i.e. each edge is a set of exactly two vertices. If an edge {x,y} e E exists, the two vertices x and y 
are connected. For our construction, we decompose the vertex set V into two disjoint subsets V\ and V2- As a 
special case, V2 may be empty. To each vertex x £ V\ we associate a complete graph K„ with n vertices, n > 2. A 
complete graph is a graph where each vertex is connected with each other vertex. K2 is an edge, Kt, is a triangle, 
K4 is a tetrahedron. These complete graphs form building blocks of the new lattice. We denote these subgraphs 
as K„{x). For a discussion of the Hubbard model on the complete graph the reader to referred to [13] and the 
references therein. 

We now construct the graph G = (V,E) as follows: the vertex set is V = V\ IJV2 where V2 — V2 and 

Vi = |J V(K n (x)) (5) 
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where V(K n (x)) is the vertex set of the complete graph K n (x) and E(K n (x)) is its edge set. The edge set E is 
constructed as follows: 

E =E X UE2UE3UE4 (6) 



£1 = |J £(£,(*)) (7) 
E 2 = {{x,y}eE:x,yeV 2 } (8) 
£3= U M:zeV(«,(i))} (9) 

x€V li yeV 2 ,{x,y}€E 

Ea= U {{z,z'}: zeV(£„M), 

xeV u x'£V u {x,x'}^E 

z' eV(K„(x'))}. (10) 

In words: if x,x* G Vi are nearest neighbours in G, we connect each vertex of K n (x) with each vertex of Knft). 
If x <G Vi, y € V2 are nearest neighbours in G, we connect each vertex of X„(x) with y. 

We associate with each edge in £, a hopping matrix element f, > and denote this new hopping matrix by T. 
We add to T on-site energies t 2 for the vertices in V 2 . The eigenstates and the spectrum of this single-particle 
hopping matrix can be constructed from the eigenstates and the spectrum of the adjacency matrix A (G) of the 
original graph G, A(G) = {a X y) xyl zv{G) where = 1 if {x,y} 6 E, otherwise. To show this, we write 

where A (y - contains the matrix elements of A(G) connecting vertices of V, to Vj. We introduce the identity 
matrices £, = (5 ZjZ /) Z)Z / 6 y. on V,- and the matrix B = (b zx ) ze v l ^eVi' ^zx = 1 if z € A"„(x) and fo^ = otherwise. B 
maps the vertices in V\ to the vertices in V\ . Then, the new hopping matrix T can be written as 



ti(BB' -E 1 )+t 4 BA n B t t 3 BA 12 

t 3 A 21 B' t 2 A 22 + t 2 E 2 



(12) 



Note that B'B = nE\, where £1 is the identity matrix on V\. We now construct all eigenstates of T. One class of 
eigenstates of T are 

W= ( o° )' ^oekernel(B') (13) 

with the eigenvalue —t\. A basis in this eigenspace is formed by all eigenstates with eigenvalue -1 of K n (x) for 
all x G Vi. K„ has one eigenvalue n — 1 belonging to the eigenstate 0q = (1,1, — , l) r and n — 1 eigenvalues -1 
belonging to eigenstates orthogonal to 0o> these are elements of the kernel of B'. The second class of eigenstates 
of T are of the form 

-(IT U )-'(£)■ 



We obtain 



^ I t3A 2 ianYi+t2A22Y2 + hY2 ' 
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where we made use of B'B = nE\. We now let a = n~ l l 2 , — h/a, ti = t^/a, and I2 = t\(n— 1). Then, we 
obtain Ty/ = [t\ [n— 1) + t4&]y/. We choose t\ and ?4 so that t\n + f4<5 m i n > where 5 m ; n is the lowest eigenvalue 
of A(G). Since 5 m j n > — d max (G), which is the maximal degree of G, it is sufficient to choose ?4 < t\n/d max (G). 
Then the ground state of T has a lowest eigenvalue — 1\ with degeneracy Nj = (n — l)|Vi| and the ground states 
are given by d 1 3t . By construction, there are |VI eigenstates with eigenvalues above — 1\ and since Nj + \V \ = \V \ 
the construction yields all eigenstates of T. 

This construction is valid even if V2 is empty. 

Let us remark that if G is a translationally-invariant lattice with r energy bands and G is translationally invariant 
as well, then G has r + n — 1 energy bands and the n — 1 additional energy bands are fiat and degenerate with 
eigenenergy — t \ . 

This construction yields a system in which the matrix p decomposes into \V\ \ blocks, one for each x€V\. The 
V x are all the same and their value is n — 1 . 

The construction can be generalized, we may choose different values of n for different x S V. 



3 Ground state properties 

We are now ready to state our main results for the Hubbard model with a lowest single-particle eigenenergy Sj 
which is A^/-fold degenerate and for which the projector onto the eigenspace of e c { fulfills the properties of Sect. 

m 

Theorem. For Hubbard models with a lowest single-particle eigenenergy e t i which is A^/-fold degenerate and 
for which the projector onto the eigenspace of fulfills the properties of Sect. [2] the following results hold for 

N e <N d : 

1. The ground state energy is E c iN e . 

2. Let A x be an arbitrary local operator, i.e. an arbitrary combination of the four creation and annihilation 
operators c\ a and c XB . The correlation function = (A x A y ) — (A x )(A y ) has a finite support for any 
fixed x and vanishes if x and y are out of different clusters V^. The system has no long-range order. 

3. The system is paramagnetic. 

4. The entropy at zero temperature S(c) is an extensive quantity, S(c) = 0(N e ). It increases as a function of 
c = N e /N c i from for c = to some maximal value S max > Y.k[( v k — l)ln2 +ln(Vjt + 2)] and then decays 
toS(l) = I ft ln(v*+l). 

Proof. We first construct a suitable basis. The different clusters Vjt are completely decoupled as long as 
ground states are considered. We therefore discuss first a single cluster. According to Q, it is possible to 
choose a single-particle basis (not orthonormal) for the space of degenerate single-particle ground states on 
the cluster with the following properties: 

1 . For each basis state y/j there exists a vertex set so that the support of y/,- is a subset of . 

2. For each basis state y>i there exists a unique x; G Vufj such that yi(jCf) > and V / i( x ;) = for all i ^ j, 
;'=!,..., \B k \. The set V k{i) \ {x h /=!,..., \B k \} is not empty. 
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Since the set Vk(j) \ {xi, i= l,...,\Bk\} is not empty, these states overlap and are not orthogonal. Let Vk = 
Using such a basis Bk, we can put electrons into the different basis states with the condition that one single- 
particle basis state contains at most one electron. In a state with two electrons in a basis state Xjfj we have a 
double occupancy on x;. Since all other basis states vanish on Xt, it is not possible to get rid of that double 
occupancy due to some linear combinations of these states. Therefore, any state with a doubly-occupied y/j has 
a non-vanishing interaction energy and cannot be a ground state. Note that the absence of doubly occupied y/j is 
a necessary but not a sufficient condition for a ground state, since it does not exclude double occupancies on the 
lattices sites in Vma \ {x;, i = 1 , . . . , \Bk | } . This construction was used in [7] to show that if there are Vk electrons 
on the cluster, they all have the same spin S% — V&/2 and that the degeneracy of the ground state on the cluster 
is 2Sk + 1 = Vk+ 1. The trivial case is one electron on the cluster, where the degeneracy is 2Vk- For electron 
numbers iik with 1 < iik < Vk it may be difficult to calculate the ground state degeneracies. But for Vk < 2 we 
have a complete description of all ground states in the cluster Vk- 

It is trivial to generalise this argument to the entire lattice using the basis Ujt^*:- The states in the different Bk 
can be filled independently. The lattices studied in [8, 9 10] all belong to the class with Vk < 2. In iflOl . it 
was assumed that the construction above yields all of the ground states and numerical results were presented to 
confirm that. Our argument is a rigorous proof of this statement. 

We now come to the four statements in the theorem. The first point is trivial. The ground state energy is e ( /N e , 
since states with that energy minimise both the kinetic energy and the interaction. The ground states have no 
doubly-occupied sites. 

For the proof of the next statements, we use a grand-canonical formulation. Let 

Z(z, xeV}) = (f° exp(£a. v A x )> (16) 

X 

be the generating function for correlation functions containing the operators A x . N e is the number operator. (•) 
denotes the ground state expectation value for arbitrary electron numbers < Nj. Since the system decomposes 
into clusters, the generating function can be written as 

z{ z , K, xev}) = Y\z(z, K, x g v k },Pk) (17) 

k 

where 

ZfcKie V k },p k ) = (f*exp(E<hAx))k (18) 

X 

denotes the ground state expectation value on the cluster V#. We have 

(Ax) = 4~ i°zfete,*e v})| ov , (19) 

aa x 

(A x Ay) = lnZ ( z ' V » U=0V* (2°) 

Since lnZ(z, {a x , x G V}) =^lnZ(z, {a x , x €Vk},Pk) one has (A x A y ) = (A x )(A y ) if x and y are out of different 
clusters Vk- Thus, Pa^cv vanishes if x and y are out of different clusters. Since this statement holds for any z, it 
holds as well for a fixed particle number N e . 

The third point follows from the second if we take for A x the local spin-operators. To be more explicit, let 
us calculate the expectation value of the total spin. It can be written as (S 2 ) = Y^ki^l) wnere $k = HxeV k ^x 
is the spin operator on the cluster k. A trivial upper bound for (Si) is lvfc(v* +2). A trivial lower bound 
for (S 2 .) is jNk, where Nk is the number of electrons on the cluster Vk. Therefore, we obtain \N e < (S 2 ) = 
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\ Lit Vk(Vk + 2) < \ Vmax(v ma x + 2)iV r . Therefore, S tol is not an extensive quantity, the system is 



StotOStot + l) < iLkVk(vk + 2) < 4 Vmax(Vi 
t ferromagnetic. The maximum value of 
2). This proves the third point in the theorem. 



not ferromagnetic. The maximum value of Stot occurs for N e = Nd, and is given by 5 to t(5t t + 1) = \ Lit v k(Vk - 



We come now to the fourth point of the theorem. We will calculate the entropy density by calculating the ground 
state degeneracy of the system. The fact that this can be done for lattice models with finite range interactions 
has been proven by Aizenman and Lieb [14]. They pointed out that there is a problem when interchanging the 
limit T — >• and the thermodynamic limit. In order to calculate the entropy density at zero temperature, one has 
to calculate the ground state degeneracy for all possible boundary conditions. This is possible in our case since 
the lattice decomposes into finite clusters, so that the boundary has no effect on the result. 

To calculate the entropy at zero temperature, let us now calculate the grand canonical partition function. Since 
the problem decomposes into a set of clusters Vk, the contribution of these multi-particle ground states to the 
grand canonical partition function is a product of the partition functions of these clusters. 



Z(z)=l\z(z,p k ). 



(21) 



The general form of Z(z,Pk) is 



z(z,p*) = £/>jv 

.7=0 



(22) 



where p v , is the number of states with j electrons on a cluster with v states. One has p$ = 1, p\ = 2v, Py = V + 1, 



'jau+i>(;) 



for 1 < j < V. The lower limit for p v - is the number of fully-polarised states with j electrons 



on a cluster with v states. 

From Z we obtain the grand canonical potential £2. The general from of £2 is 



a = -j8- 1 ^lnZ(z,p (fe ). 



(23) 



The entropy is 



S(z) = ~ = £lnZ( Z ,p*) - z lnz£^lnZ( Z ,p fc ). 



' dz 



(24) 



k k 

Let us introduce c(z) = N e (z)/N r . Since N e (z) = ~£ we obtain c(z) = ^ £k N k(z) where 

Nkiz) = Z(z,p k ) dz 

is the number of particles on the cluster VV One has S(z) — Y,kSk(z), Sk(z) = lnZ(z,pk) — A^(z)lnz. 
c{z) is a strongly monotonically increasing function of z- For the derivative of the entropy, we obtain 



(25) 



dS 

dz 



= — lnz 
= -lnz 



^ d\n.Z{z,Pk) ^d 2 lnZ(z,pk) 
£ dz 



dz 1 



1+ — 

dz 



-Wz) 



(26) 



-AT, lnz 



dc 
~dz 
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Since 37 > 0, the only maximum of 5(z) occurs m at z = 1. The value is S(z = 1) = J^lnZ^pjfc). F° r small 
values of z we have 

Z(z,p k ) = l+2v k z + 0(z 2 ) (27) 

and therefore 

S(z)=2z(l-lnz)N d (l+0(z)) (28) 

In the limit z — >• °°, at the maximal density, the degeneracy in the cluster p,t is + 1 . The total degeneracy is 
rifc( v * + 1) anc l tne entropy is S(z — > <*>) = L<.-l n ( v A. + 1). Therefore, the entropy increases monotonically to it's 
maximum at z — 1, 5(z = 1) = £j.lnZ(l,p4;) and then decays monotonically to 5(z — > °°) = L^ln(v^ + 1). Using 

pj > (7 + 1) ^ ^ ^ we obtain the lower bound for 5(z = 1) in point 4 of the theorem. Since c(z) is strongly 

monotonically increasing, these properties hold for 5(c) as well. 

3.1 Examples 

For clusters with = 1 one has 



Z(z,p k ) = l+2z. (29) 



For V£ = 2 one obtains 



ZfcpA-) = l+4z + 3z 2 (30) 



since there is one state on the cluster Vk with no particles, four states with one particle and three states with two 
particles. For Vk = 3 one has 2v k = 6 states with one electron and Vk + 1 = 4 states with three electrons. Two 
electrons on the cluster can form a triplet state or a singlet state. One gets 3v k —9 triplet states and between 
and 2 singlet states. In the basis Bk none of the basis states can be doubly occupied. Therefore one has only 
three different pairs which could form a singlet. But, since pk is irreducible, at most two different pairs without 
a doubly occupied site can be constructed. Therefore, for v k = 3 we obtain 

Z(z,Pk) = l+6z+(9 + s k )z 2 +4z 3 (31) 

where s k is the number of possible singlets on the cluster p k , which can be 0, 1, or 2. 
If V k < 3 for all k, the total partition function is thus 

Z(z) = (l+2z) N '-°(l+4z + 3z 2 f 2 ' 
(l+6z + 9z 2 + 4zy 3 - 

(l+6z+10z + 4zY 3 'i (32) 
(l+6z+llz 2 + 4z 3 ) A '" 

where N v , s is the number of clusters with Vk = V and s possible singlet states. From d32l one obtains the grand 
canonical potential 

^ = -j3- 1 (« 1 , ln(l+2z)+n 2 ,oln(l+4z + 3z 2 ) 



+n 3)0 ln(l +6z + 9z 2 + 4z 3 ) 

+« 3 ,iln(l+6z+10z 2 + 4z 3 ) (33) 



+» 3 ,2ln(l+6z+llz 2 + 4z 3 )) 
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where n VJ —N v , s /N r , «i,o + «2,o + »3,o + »3,i + «3.2 = 1- 

Let us mention that all of the one-dimensional lattices treated in iflOl belong to this class with «2,o = 1, i.e. 
Ni i0 = and N 3iS = for s = 0, 1,2. 

Inverting c(z) we obtain z as a function of c and then S(c). There are two cases where c(z) can be inverted 
easily, the case where all Vj. = 1 and the case where all — 2. In the first case one obtains a linear equation 
for z(c), in the second case a quadratic. The one-dimensional models in [10] are a special case for — 2. The 
result given there is valid for any model where all clusters have = 2. For the case < 3, ( 1331 ). one obtains 

2z 4z + 6z 2 

c = "1,0 r , — + »2,0 
+«3,0 



l+2z /AJ l+4z + 3z 2 
6z+18z 2 +12z 3 



l+6z + 9z 2 + 4z 3 
6z + 20 z 2 +12z 3 
10z 2 + 4; 
6z + 22z 2 +12z 3 



+ ^ 1 l+6z+10z 2 + 4z 3 (34) 

+"3,2 



l+6z+llz 2 + 4z 



One has c < n\ + 2«2 + 3(«3,o + »3.i +"3.2) = 1 + n 2 + 2(«3.o + «3,i +13 2), where the upper limit is reached for 
z — >• 00. In this limit, the entropy takes the value 



5 =A^ r («i.o + 2« 2 ,o + 3(n3.o + n3.i +« 3 , 2 )). (35) 



The maximum 



of the entropy occurs at 



5 = A^ r (ni,oln3 +«2,oln8 +«3 : oln20 

+ n 3 ,iln21+«3 )2 ln22) (36) 

2 5 9 38 5 

C = -Ml,0 + ^"2,0 + ^«3.0 + 2]" n 3,l + 4»3,2- (37) 

For the examples constructed in Sect. [2] the construction of this basis is easy. We simply use the states with the 
properties = 1, y/,(«) = — 1, YiU) = for j =^ i,n, i = 1, ... ,n — 1. The partition function is 

Z(z) = J] Z(z,K„(x))=Y\Z(z,K n ) N ". (38) 

It is sufficient to consider a single K„. Let V(X„) = {1 ,...,«}. The cases n < 3 correspond to the cases < 2 
already discussed above. ^4 has = 3. For it is possible to construct two pairs of non-overlapping single- 
particle ground states. Using the basis introduced above the two pairs are y/2 — 1/3} and {l/Aj — y/2, 1^3}- 
Therefore, we have A^o = ^31 = in that case, and Z(z,Kn) = 1 +6z+ 1 lz 2 +4z 3 . For larger values of n the 
number of non-trivial cases increases rapidly. In principle it is possible to completely describe the multi -particle 
states as well. We do not discuss these cases here. 

The most important point in the discussion of these systems and of the multi-particle ground states is that 
although this construction allows for many different examples of solvable systems in arbitrary dimensions (since 
in the above construction G may be an arbitrary lattice in arbitrary dimensions), the ground state properties for 
T = and N e < and the contribution of the ground states to the low temperature properties of the system are 
that of a collection of zero-dimensional systems. The properties do not depend on the dimension of the lattice, 
but only on the number of different subgraphs of type K n the lattice contains. 
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4 Summary and Outlook 



This paper yields a complete classification of all Hubbard models for which the degeneracy Nj of the single- 
particle ground states is some finite fraction of the number of lattice sites and for which the projector onto this 
subspace is highly reducible, i.e. where the number of irreducible submatrices N r of this projector is some 
finite fraction of the number of lattice sites. Each subspace lives on a local cluster and different clusters do not 
overlap. We show how lattices with these properties can be constructed in arbitrary dimensions and we derive 
some properties of the ground states of such models for electron numbers N e < N4. Examples of such lattices 
in one and two dimensions were previously presented by Batista and Shastry (8), Maksymenko et al. [10|, 
and others, see also the references therein. Maksymenko et al. [10] gave a rather complete discussion of the 
Hubbard model on some one-dimensional lattices of this type. 

The important point is that the ground states properties for such models do not depend on the details of the 
lattice or on its dimensionality, but only on the properties of the local clusters. Global properties like the 
entropy at T = can be calculated. The behaviour of the entropy is similar for all of these lattices. The entropy 
density as a function of the density of particles grows from to some maximum and then decays to some 
finite value at N e — N c /. Thus, the one-dimensional lattices in Jl0| are ideal prototypes of all these models, 
and no essentially new physics occurs in the higher-dimensional models. Maksymenko et al. [ 1 1 discussed 
the case where the degeneracy within the cluster is lifted. In that case the system still has a large ground 
state degeneracy, finite entropy density, etc. In that case, the ground states still are located on the small local 
clusters and the dimensionality of the lattice remains unimportant. This may of course change if one lifts the 
degeneracy by some small perturbation so that the lowest bands are no longer strictly flat. For the discussion 
of the stability of ferromagnetism the situation then becomes much more difficult, see HE]. We expect that 
with such perturbations, the dimensionality of the lattice becomes important as well and that the analysis will 
be much more difficult. Nevertheless, in that case new and interesting physics may occur. 

Maksymenko et al. \ 1 1 discussed not only the Hubbard model but also the Heisenberg model on such lattices. 
We expect that their results for the two one-dimensional models can easily be generalised to the class of lattices 
described here as well. We expect that, as for the Hubbard model, one obtains no new interesting physics. 

Another class of models which are closely related to spin systems with antiferromagnetic exchange interactions 
are bosonic Hubbard models with flat bands on similar lattices, see e.g. ifTTl [TBI and the references therein. 
Bosonic Hubbard models can also be discussed on the lattices presented here with similar results. 

In this paper, we only discussed the T — properties of these models. For T > 0, the situation becomes more 
complicated. The detailed structure of the lattice and the dimensionality become important, since the other 
single-particle eigenstates ( TT4l > depend on the detailed lattice properties. 
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